Provenance and family variations in early growth of Manchurian walnut (Juglans mandshurica Maxim.) and selection of superior families

This study, conducted in China in November 2020, was aimed at exploring the variations in growth traits among different provenances and families as well as to select elite materials of Juglans mandshurica. Thus, seeds of 44 families from six J. mandshurica provenances in Heilongjiang and Jilin provinces were sown in the nursery and then transplanted out in the field. At the age of 5 years, seven growth traits were assessed, and a comprehensive analysis was conducted as well as selection of provenance and families. Analysis of variance revealed statistically significant (P < 0.01) differences in seven growth traits among different provenances and families, thereby justifying the pursuit of further breeding endeavors. The genetic coefficient of variation (GCV) for all traits ranged from 5.44% (branch angle) to 21.95% (tree height) whereas the phenotypic coefficient of variation (PCV) ranged from 13.74% (tapering) to 38.50% (branch number per node), indicating considerable variability across the traits. Further, all the studied traits except stem straightness degree, branch angle and branch number per node, showed high heritability (Tree height, ground diameter, mean crown width and tapering, over 0.7±0.073), indicating that the variation in these traits is primarily driven by genetic factors. Correlation analysis revealed a strong positive correlation (r > 0.8) between tree height and ground diameter (r = 0.86), tree height and mean crown width (r = 0.82), and ground diameter and mean crown width (r = 0.83). This suggests that these relationships can be employed for more precise predictions of the growth and morphological characteristics of trees, as well as the selection of superior materials. There was a strong correlation between temperature factors and growth traits. Based on the comprehensive scores in this study, Sanchazi was selected as elite provenance. Using the top-percentile selection criteria, SC1, SC8, DJC15, and DQ18 were selected as elite families. These selected families exhibit genetic gains of over 10% in tree height, ground diameter and mean crown width, signifying their significant potential in forestry for enhancing timber production and reducing production cycles, thereby contributing to sustainable forest management. In this study, the growth traits of J. mandshurica were found to exhibit stable variation, and there were correlations between these traits. The selected elite provenance and families of J. mandshurica showed faster growth, which is advantageous for the subsequent breeding and promotion of improved J. mandshurica varieties.


Introduction
The genus Juglans is one of the most important genera of Juglandaceae family, comprises about 21 species, and Juglans mandshurica Maxim. is a deciduous fruit tree belonging to the Juglandaceae family [1,2].This tree is native to northern and northeastern China and is mainly distributed in the Changbai and Xiaoxingan Mountains at an altitude ranging from 500 m to 1000 m [3].The species also grows naturally in Russia, Japan, and Korea and can tolerate temperatures of -50˚C, and has important economic values [4,5].The J. mandshurica is a versatile wood that is rigid, wear and corrosion-resistant, easily processed, and produces high quality wood suitable for joinery in luxurious furnishings [6,7].Besides, its seeds are rich in nutrition and are considered a premium forest food product [8].Due to its high economic returns and nutritional value, its production and consumption are increasing all over the world [9,10].Furthermore, J. mandshurica has high medicinal and health values [11].For instance, its husk, bark, roots and leaves contain juglone, which has good antitumor [12], antifungal [13], antiviral [14], antioxidant [15], anthelmintic [16] and hypoglycemic activities [17].They are widely used in pharmaceutical, food, and cosmetic industries.
Understanding and leveraging the variation and interplay of genetic and phenotypic traits hold decisive importance in optimizing the improvement of biological characteristics across various research domains.Against this backdrop, techniques such as correlation analysis, general combining ability analysis, phenotypic coefficient of variation analysis, genetic coefficient of variation analysis, and family heritability analysis have been widely employed to interpret and exploit genetic diversity.Prior studies have laid the foundation for plant breeding by exploring genetic associations between traits, evaluating parental combining effects, and analyzing phenotypic and genetic variations [18,19].Recently, researchers have delved deeper into understanding the specific roles of genotype and environmental factors in phenotypic expression and genetic variation using more advanced and refined methods [20].While these studies have offered valuable insights, there remain unresolved issues and challenges in the specific application of Juglans mandshurica.
Provenance tests are crucial for species that have high economic or ecological value [21].Systematic provenance test can explore the adaptability of different provenance materials to the local environment, quantify the genetic population differentiation, and provide materials for forest breeding and ex-situ protection [22,23].Provenance test for J. mandshurica began in the 1980s ~1990s [24] by investigating and analyzing phenotypic growth traits of eight J. mandshurica provenances, which resulted in selection of two best provenances in Kuandian and Shulan.Based on a geographical and climatic variation study, Liu profiled the differences in growth and physiological performance of the provenances, and preliminarily classified J. mandshurica into four major provenance zones [25].Since then, there have been several studies on the provenance of J. mandshurica in China [8,[26][27][28][29][30][31][32][33][34].Most of these studies, however, concentrated on seed traits or a few growth traits, other growth data of J. mandshurica seedlings are relatively scarce.In addition, due to challenges posed by late genetic improvement of the J. mandshurica, coupled with long sexual reproduction cycle and low survival rate of vegetative propagation such as cuttings and grafting.These led to a small number and usage rate of improved varieties of J. mandshurica.In addition, excessive logging and man-made destruction of natural forests in recent decades have led to a large reduction in the number of trees [32].Some available natural distribution areas are even mostly secondary forests, and J. mandshurica trees are on the verge of depletion [35].At the same time, climate and environmental changes include warming temperatures, late spring frost, drought stress, high temperature stress and salt stress also negatively affected the growth of genus Juglans [36][37][38].It is, therefore, important to explore variation patterns of provenances and families to guide collection and conservation of J. mandshurica resources, protection of biodiversity and maintenance of ecological balance.Recent finding shown that although J. mandshurica natural forest is seriously damaged, there is still a high degree of genetic diversity [6].The result showed that the existing resources of J. mandshurica may still have rich variation, which provided the possibility for the selection of improved varieties of J. mandshurica.On the other hand, the main breeding goal of genus Juglans plants is to cultivate productive new varieties and rootstocks with superior resistance [39][40][41][42][43][44].Due to its good drought and cold tolerance and well-developed root system, J. mandshurica is often used as an excellent rootstock for walnuts, which can effectively improve walnut yield and stress resistance [45,46].Therefore, preliminary selecting of J. mandshurica resources can provide basic materials for the breeding of high-yield and high-resistance J. mandshurica varieties.This study employed the provenance test to profile 5-year-old J. mandshurica tress in Wanrenhuan forest farm of Binxian county, Heilongjiang Province.The objectives were to: (1) evaluate the variations among different J. mandshurica provenances and families.(2) select elite provenance and families.The existing studies have primarily focused on preliminary insights into the economic, medicinal, and timber value of J. mandshurica, as well as the presence of certain bioactive components.However, these studies have limitations in the in-depth exploration of diversity, the selection of genetic sources, and the evaluation of families for J. mandshurica.They also do not provide comprehensive information on growth characteristics.Therefore, the primary strength of this study lies in the indepth investigation of the growth traits of J. mandshurica.It involves the assessment of genetic variations and the selection of elite provenance and families, thereby providing robust support for the forest management and timber production of this species.

Experimental site and plant materials
Seeds of 44 half-sib families (44 individual trees, each tree collected 120 seeds) from six J. mandshurica provenances in Heilongjiang and Jilin provinces, China, were collected.The provenances were Sanchazi (15 families), Hulin (3 families), Daquanzi (9 families), Dongjingcheng (10 families), Yabuli (4 families) and Tieli (3 families) (Table 1).The seeds from open-pollination were sown in autumn 2015, with a sowing area of 1.69 hm 2 , each family selects 100 healthy seeds for sowing in the nursery, with a germination rate of 74%.32 seedlings with consistent growth of J. mandshurica were selected per family for afforestation, planted the following year in a randomized complete block design, with 8 trees and 4 blocks at a spacing of 3 m × 4 m (plant spacing and row spacing) in Wanrenhuan Forest Farm of Binxian County, Heilongjiang province (E45˚44 0 20@ ~45˚57 0 18@, E128˚04 0 29@ ~128˚19 0 13@), and the distance between families is 4 m.The altitude of experimental area is 340 m, with a mean annual temperature of 2.4˚C, mean annual rainfall of 577 mm, mean annual sunshine hours of 2090 h, and a frostfree period of 130 d.The geographical location and environmental factors of the provenances are given in Fig 1

Trait measurements
In November 2020, we assessed the growth and stem traits of a total of 1408 individual tree (The test forest consists of 44 families, with 32 individual trees per families) in a 5-year-old J. mandshurica provenance test forest.The traits included tree height (m), ground diameter (cm), mean crown width (m), stem straightness degree, tapering, branch angle (˚) and branch number per node.The tree height was measured using a Box staff.The ground diameter was measured using a tree caliper.The east-west and north-south crown width were measured using a stainless-steel measuring tape, and the average of two crown width was taken as the mean crown width value.The stem straightness degree was meas ured using grading method and a score from 1 to 5 was assigned, where score 1 was assigned if the stem had more than two obviously bend points in the stem and a score of 5 if the stem was completely straight [47].Stem straightness degree necessitates the square root transformation of data during the process of conducting analysis of variance (ANOVA).The ratio of ground diameter to tree height was used to calculate tapering [48].The branch angle was measured using a protractor.Climatic factors were obtained from the National Meteorological Information Center (http://data.cma.cn/).In this study, field site access and sample collection did not require specific permits.The research was conducted with the full consent and cooperation of the relevant site managers and authorities.No permits were necessary due to the non-sensitive nature of the work and the site's accessibility to researchers.In addition, we collected growth measurement data from the Experimental Forest.These data, referred to as 'Experimental Forest Growth Measurement Data,' have been included as supporting information in this study.

Statistical analysis
Mixed linear model was applied to investigate into the differences of provenance, family and block effects on growth traits.The linear model was applied using the following formula [49].All data were tested for normal distribution and homogeneity of variances.where X ijkl is the observed value of an individual tree l in family j within provenance i growing in block k, μ is the family mean, B k is the fixed effect of block k, P i is the fixed effect of provenance i, F j(i) is the random effect of family j within provenance i, BP ik is the fixed interactive effect of block k and provenance i, BF j(i)k is the random interactive effect of block k and family j within provenance i, and e ijkl is the random error.Heritability (h 2 ) among families for each trait was calculated following a formula described by Zhang [50]: where s 2 F ; s 2 FB and s 2 e denotes the variance components of family, family by block interaction, and residual error, respectively.B and N represent the number of blocks and replicates, respectively.
The coefficient of variation was calculated by using the following formula [51]: phenotypic coefficient of variation (PCV): genetic coefficient of variation (GCV): where s 2 P denotes the phenotypic variance of a trait, s 2 g denotes the additive genetic variance component in maternal half-sib of a trait (4s 2 F ), � X is the mean value of a growth trait.The Pearson correlation coefficient was used to measure the linear relationship between two variables.The Pearson correlation coefficient (r) is a standard statistical method employed to quantify the correlation between two variables.Its calculation formula is as follows: r AðxyÞ ¼ s aðxyÞ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi s 2 aðxÞ � s 2 where σ a(xy) is the covariance component between trait or environmental factor x and y, s 2 aðxÞ and s 2 aðyÞ are the variance component for trait or environmental factor x and y, respectively.The general combining ability (GCA) was calculated by the formula [52]: where g is the general combining ability of parents, x is the offspring mean of a hybrid combination of the parents for a particular trait, and μ is the offspring mean of all combinations for this trait.
The principal component values and comprehensive score were calculated according to Chen [52]: where Y i is the value of principal component i, α ij is the eigenvalue of trait j within principal component i, X j is the mean of trait j, W is the comprehensive score, ω i the contribution rate of principal component i, n is the number of traits, and p is the number of extracted principal components.
Selecting elite provenance and families with a selection rate of 10%.The real gains of elite provenance and genetic gain of elite families were computed using the following formula [49]: genetic gains (ΔG 1 ): real gains (ΔG 2 ): where h 2 is the heritability of family, W the selection differential, � X is the family or provenance mean.
The mean analysis, ANOVA, principal component analysis and Pearson correlation analysis were conducted using IBM SPSS Statistics, version 26.0 (IBM Corp., Armonk, NY, USA) [53].Additionally, variance parameters and general combining ability were calculated.For all statistical analyses, P-values less than 0.05 were considered statistically significant.Provenance distribution map was drawn using ArcGIS, version 10.7 (Esri, Redlands, CA, USA) [54].Custom data visualization was carried out in the R environment, version 4.0.2, using "corrplot" packages (R Core Team, Vienna, Austria) [55].

Variations in early growth among provenances and families
Significant variations (P<0.01) in growth traits were detected among different provenances and families within provenances (Table 2).Description statistics of growth and stem traits together with GCV, PCV and heritability are shown in Table 3.The GCVs of all the traits ranged from 5.44 to 21.95%, among which tree height, ground diameter and mean crown width were higher than 15%.The PCVs of all the traits ranged from 13.74% to 38.50%, among which the stem straightness degree and branch number per node were more than 30%.The heritability of all traits in the families ranged from 0.260 to 0.908, among which the heritability of all the traits was higher than 0.7, except for the stem straightness degree, branch angle and branch number per node.In addition, the mean value of all growth traits for different provenances and families are shown in S1 and S2 Tables.The results in S1 Table showed that Daquanzi and Sanchazi provenance had more elite growth traits, which was also reflected in S2 Table .Most of the families with elite growth traits from the Daquanzi and Sanchazi provenance.

Correlation analysis
The correlation among the tree height, ground diameter and mean crown width were significant and positively correlated with a coefficient exceeding 0.8 (Fig 3).Stem straightness degree, branch angle and branch number per node had a weak positive correlation with tree height, ground diameter and mean crown width, with a correlation coefficient of not more than 0.3.Among environmental factors, annual temperature had a strong positive correlation with mean crown width, while annual temperature had a negative correlation with tapering, with correlation coefficients of more than 0.8 (Fig 4).Additionally, it was observed that January mean temperature and winter mean temperature had a strong positive correlation with tree height, ground diameter and mean crown width, the magnitude of the correlation coefficients exceeded 0.9.Tapering was strongly positively correlated with altitude (r = 0.84).In addition, stem straightness degree was strongly positively correlated with winter mean rainfall (r = 0.84).No correlation was observed between the remaining traits and the environmental factors.

Principal components analysis
We performed principal components analysis (PCA) and three principal components were extracted based on eigenvalues greater than 1 (Table 5 and Fig 5).The eigenvalue for PC1 was 2.826 and the contribution rate was 40.37%.The absolute eigenvalues for the tree height, ground diameter and mean crown width were relatively high at 0.934, 0.942 and 0.920, respectively.On the other hand, the eigenvalue for PC2 was 1.154 and the contribution rate was 16.49%.The absolute eigenvalue for tapering was relatively high, at 0.926.The eigenvalue for PC3 was 1.007 while the contribution rate was 14.39%.We found that the absolute eigenvalues  for the branch angle and branch number per node were relatively high, at 0.623 and 0.578, respectively.The cumulative contribution rate of the three principal components was 71.25%, which covered most of the information of growth and stem form traits for the tested families.

Selection of elite provenance and family
The linear equation for each principal component score was obtained from the results of principal component analysis, and Y 3 ¼ À 0:042x 1 À 0:163x 2 À 0:059x 3 À 0:436x 4 À 0:251x 5 þ 0:623x 6 þ 0:578x 7 , where x 1 , x 2 , x 3 , x 4 , x 5 , x 6 and x 7 denote the tree height, ground diameter, mean crown width, stem  straightness degree, tapering, branch angle and branch number per node, respectively.The contribution rate of each principal component was used as the weight to calculate the comprehensive score, using the following formula: The final comprehensive scores and ranking of provenances and families are given in Tables 6 and  7.The provenance was selected based of tree height, ground diameter and mean crown width, with a real gain of 5% or more.The real gain for each trait was 6.46%, 6.73%, 7.31%, 4.11%, 0.04%, 0.01% and 0.47% for tree height, ground diameter, mean crown width, stem straightness degree, tapering, branch angle and branch number per node, respectively (Fig 6).We then applied 10% inclusion rate and selected four families: SC1, SC8, DJC15 and DQ18.The mean values of the selected elite families for tree height, ground diameter, mean crown width, stem straightness, tapering, branch angle, and branch number per node were 2.31 m, 5.91 cm, 1.60 m, 4.04, 2.58, 49.85˚and 2.20, respectively.The genetic gain achieved after selecting these superior families were 10.47%, 13.51%, 14.83%, 3.51%, 0.94%, 5.49% and 1.72%, for tree height, ground diameter, mean crown width, stem straightness, tapering, branch angle, and branch number per node, respectively (Fig 6).

Discussion
Genetic variation drives evolution of different species, thereby enabling them to adapt to the prevailing environmental conditions [56].Thus, assessment of genetic variation in different populations is important as a first step in breeding improved genotypes and for effective conservation of forest genetic resources, and the results will be helpful for the breeding of highquality commercial varieties [57,58].The variation in phenotypic traits is the result of the interaction between genetic variation and environmental factors, reflecting the variability of genotypes, populations, ecotypes, and provenance [59].Our results demonstrated that most of the traits varied significantly among different provenances and families, while the differences among blocks were small, which was consistent with data from previous studies [27,33,34,60,61].Similar research results were also obtained in Eucalyptus urophylla and some subtropical pine species [62,63].The results also demonstrated that the accessions had comparatively large variation for genetic improvement, and the differences mainly resulted from genetic rather than environmental stimuli.It should be noted that in natural environments different populations often have significant differences due to interaction between genetic and environmental diversity.Thus, analyses of the variations have contributed to a better understanding of the inherent patterns of genetic variation [64].
In this study, the mean values of tree height and ground diameter for 5-year-old trees were higher than previous findings that used 7-year-old J. mandshurica [33] and 5-year-old J. mandshurica [65,66].The accessions showed a distinct growth advantage.In comparison to the 1-9-year-old walnut trees in the edge habitats of natural forests, our research material did not exhibit a significant advantage in tree height growth [67].This phenomenon appears to be influenced by a variety of interrelated factors.Notably, the high-density planting method in artificial forests may exacerbate competition among trees, thereby imposing certain limitations on their growth.Additionally, it is worth noting that nutrient supply in artificial forests is typically relatively inadequate, necessitating supplementary fertilization and management interventions.Moreover, adverse growth conditions such as limited light exposure and water availability should be considered, as they may significantly impact tree height.By enhancing management practices and nutrient supply, there is the potential to elevate the tree height of walnut trees in artificial forests, gradually reducing the disparity in tree height compared to natural forests.The crown is a major primary site of photosynthesis and energy retention for trees, which reflects the viability and competitiveness of tree growth [68].The mean crown width in this study were 1.33 m, which was slightly lower than the results by Zhao on a 6-yearold J. mandshurica [69].This variation could be associated with the variety used and the soil fertility [70,71].In addition, there are few studies on the traits reflecting the stem form characteristics such as stem straightness degree, tapering, branch angle and branch number per node of the J. mandshurica.Therefore, our study investigated and analyzed these traits to gain a deeper understanding of the growth variation pattern of the J. mandshurica.The form quality of the trees is controlled by their own genetic and environment conditions, which could directly affect the processing, quality, yield, and commercial value of the wood [72].Larger stem straightness degree and smaller tapering have been shown to be beneficial for the selection of large-size timber, while larger branch angle and branch number per node could be conducive in improving fruit production [73,74].ANOVA results showed that the significant differences in stem form traits among the different provenances and families in this study, which could provide materials for directional breeding of J. mandshurica for diverse objectives.
The coefficient of variation is a small-scale variability measure used to reflect the genetic variability of traits in population, and the larger the coefficient of variation, the more favorable the selection of the elite materials [75,76].In this study, the phenotypic coefficient of variation (PCV) and genetic coefficient of variation (GCV) were higher than a previous report by You in Hebei province [30], thus indicating obvious differences in geographical provenance.Among the analyzed traits, the coefficient of variation was relatively large for tree height, ground diameter, mean crown width and stem straightness degree, demonstrating that using these traits as the selection index could yield more selection potential.In addition, the genetic coefficients of variation for traits such as tree height, ground diameter, mean crown width, and tapering accounted for more than 65% of the phenotypic coefficient of variation, which was consistent with the findings for Pinus koraiensis [77].This data indicated that growth variation among families is more controlled by genetic factors, which provides a basis for the selection of elite provenances and families.Besides, heritability is an important indicator in forest breeding research, which reflects the relative roles of genetics and environment in the expression of various traits, while helping to rank the importance of each trait in hybrid breeding, and provide support for early selection of trees [78,79].In this study, except for the stem straightness degree, branch angle and branch number per node, the heritability of all the traits in the family was high and could be inherited more stably.The outcomes were better than for 4-year-old, 5-year-old, 6-year-old and 7-year-old J. mandshurica [61,65,80].It has been shown that high genetic control could maximize the genetic gain and facilitate the selection of elite families [81].
In tree breeding, it is often desired to improve comprehensive traits in the target material [82].Thus, the correlation between the traits is particularly important [83].Correlation analyses and coefficients could effectively and quantitatively characterize the degree of association among traits and provide space for selection of comprehensive traits [84].However, data on the correlations among traits in the provenances of J. mandshurica remains limited.Here, we analyzed the correlation among and between traits as well as environmental factors in juvenile J. mandshurica.The data showed large correlation coefficients (r > 0.8) among some growth traits.This finding was similar to that observed in a previous study with growth traits of P. sylvestris clones [85], which indicated strong associations among these traits and were favorable for comprehensive selection.The correlation coefficients among stem straightness degree, tapering, branch angle and branch number per node were small (r < 0.3), and weakly correlated, as in the case of previous findings on P. koraiensis and P. Sylvestris [77,86].These data demonstrate that the genetic independence between growth and stem form traits and could be used for the evaluation and selection based on different breeding objectives to improve breeding efficiency.Furthermore, it was also observed in this study that the January mean temperature and winter mean temperature has a strong positive correlation with the growth traits.This indicates that the temperature in cold season has a significant effect on the growth of J. mandshurica at seedling.Similar results were reported by Liu, who studied P. sylvestris.The investigation may provide references for cultivation and promotion of J. mandshurica [87].
Phenotypic variation is an effect of genetic and environmental interactions, and environmental factors play a crucial role in shaping the plant phenotypes, such as phenotypes and nutritional traits in the fruit [88,89].The correlations between most traits and geographic factors, such as longitude, latitude and altitude, in our study were insignificant, which is consistent with results by Xia [29].The study demonstrated that the growth of different J. mandshurica provenances is weakly associated with geographic factors and their changes show a random variation.On the other hand, annual temperature had a strong positive correlation (r > 0.8) with mean crown width.This finding was consistent with the findings on Xanthoceras sorbifolium [90], which indicated that the increase in temperature is conducive to physiological activities such as photosynthesis of J. mandshurica, thus promoting the growth and development of the crown.In addition, the annual temperature showed a strong negative correlation (r = -0.85)with tapering, indicating that the higher the temperature, the relatively more use of biomass, and the tree height growth and tree morphology tends to be slender and tall, which is different from studies that evaluated coniferous tree species such as Larix Mill.[91] and might be due to the difference of the tree species.
Combining ability is one of the breeding objectives, which refers to the relative potential of parental dominance hybrid to progeny [92].It is widely used in the breeding of cross-pollinated plants, and the data could help in the selection of elite hybrid parents [93].In this study, the data showed that the general combining ability of different traits varied widely among different families, which was similar with the results from P. koraiensis [94], it was difficult to perform comprehensive selection.Therefore, it was necessary to carry out further selection in combination with PCA.PCA is a common multivariate analysis tool that helps reduce data dimensionality and retain data trends and patterns [95], and is widely used in trait analysis of walnut tree species [96][97][98].In our study, the higher eigenvalues for tree height, ground diameter, mean crown width in PC1 could reflect growth patterns.In addition, the higher eigenvalues for tapering in PC2 could represent stem form, while that of branch angle and branch number per node in PC3 could represent the branching features.This result is consistent with the data from Yang on Castanopsis hystrix [99], suggesting that the differences in J. mandshurica family are mainly in growth, stem form and branching.Based on the PCA analysis, we primarily selected one elite provenance (Sanchazi) and four elite families (SC1, SC68, DJC15 and DQ18).The selected families, whose genetic gain in tree height was higher than the results of Chu which used 15-year-old J. mandshurica [31], had obvious advantages in growth and had high genetic gain.The differentiation of materials in our study was large, which has a large potential for genetic improvement and is conducive for early selection new J. mandshurica varieties.The selected families could provide material for regional afforestation.

Conclusions
Taken together, our study successfully investigated the variations in growth and stem form traits of 5-year-old J. mandshurica provenances and families within provenances, revealing significant genetic factors at play.We identified elite families, with particular focus on tree height, ground diameter, and mean crown width.This provides a robust foundation for comprehensive multitrait selection, offering potential for the breeding and promotion of J. mandshurica varieties.This research contributes to the theoretical and breeding basis for the conservation and utilization of valuable resources, as well as the development of the local J. mandshurica cultivation industry.